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Non-isothermal fluctuating hydrodynamics and Brownian motion 
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Institut fur Theoretische Physik, Universitdt Leipzig, Postfach 100 920, D-04009 Leipzig, Germany 

The classical theory of Brownian dynamics follows from coarse-graining the underlying linearized 
fluctuating hydrodynamics of the solvent. We extend this procedure to globally non-isothermal con¬ 
ditions, requiring only a local thermal equilibration of the solvent. Starting from the conservation 
laws, we establish the stochastic equations of motion for the fluid momentum fluctuations in the 
presence of a suspended Brownian particle. These are then contracted to the non-isothermal gen¬ 
eralized Langevin description of the suspended particle alone, for which the coupling to stochastic 
temperature fluctuations is found to be negligible under typical experimental conditions. 

PACS numbers: 05.40.Jc, 05.70.Ln 


I. INTRODUCTION 

The microscopic equations of motion for strongly in¬ 
teracting many-body systems are, in general, intractable. 
A notable exception is provided by systems exhibiting 
a scale separation that allows for major simplifications 
of these equations, making them practically (and not 
only formally) useful in a wide range of experimental 
and technological applications. Of particular relevance 
is the so-called Brownian motion of a reduced number 
of slow degrees of freedom, for which the many fast de¬ 
grees of freedom provide an effective background noise. 
As Einstein realized early-on [l|, the crucial simplifica¬ 
tion arises from the assumption that the microscopic and 
Brownian degrees of freedom are in thermal equilibrium, 
which allows for a universal characterization of the noise 
dynamics without explicit microscopic calculations. The 
corresponding theory of isothermal Brownian motion is 
by now firmly established and usually additionally ex¬ 
ploits the fact that the mesoscopic degrees of freedom 
mediating between the Brownian scale and the micro¬ 
scopic noise degrees of freedom admit a coarse-grained, 
hydrodynamic description, without loss of generality. In 
particular, starting with early work by Zwanzig Q, sev¬ 
eral papers have explicitly derived (generalized) Langevin 
equations describing Brownian motion as a contraction 
of the more detailed description of a fluid governed by 
linear fluctuating hydrodynamics M. Among the ma¬ 
jor outcomes of this inquiry there is the explanation of 
the long-time tails in the Brownian velocity autocorrela¬ 
tion function Q and the robustness of the fluctuation- 
dissipation theorem against variations of microscopic de¬ 
tails and even hydrodynamic specifications, such as the 
(in)compressibility of the solvent [1, 0 or its (no- )slip 
boundary condition at the Brownian particle surface [lOj . 
The necessary nanotechnological tools to conduct quan¬ 
titative experimental tests of these ground-breaking the¬ 
oretical developments have only become available very 
recently |lll-ll3l| and vindicated the central theoretical 
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premise —i.e. the assumption of an underlying isother¬ 
mal fluctuating solvent hydrodynamics— with impressive 
precision. 

Conversely, very little is known about Brownian mo¬ 
tion in non-equilibrium solvents where the validity of a 
Langevin description is not a priori ensured and stan¬ 
dard recipes to leapfrog the microscopic dynamics using 
results from equilibrium statistical mechanics, such as 
energy equipartition, are not available. Yet, micro- and 
nanoscale motion under non-equilibrium (and in partic¬ 
ular non-isothermal) conditions is becoming increasingly 
relevant for innovative experimental and nanotechnolog¬ 
ical applications [l3 - [l^ . 

Linear fluctuating hydrodynamics, originally intro¬ 
duced by Landau and Lifshitz to describe density, mo¬ 
mentum and energy fluctuations of a fluid in a global 
equilibrium state , was later extended to non¬ 
equilibrium conditions, e.g. when a temperature gradient 
is present [isl - l^ . The efficacy of this non-equilibrium 
theory in describing fluid fluctuations is testified by the 
equivalence of its predictions to those of kinetic theory 
[2^ —within its range of validity, i.e. for dilute gases— 
and mode-coupling theory [234^ . and by the good 
agreement with light-scattering experiments; see Ref. 
[25l| for a review. In view of this success, one may ex¬ 
pect the theory to be as effective in deriving reduced de¬ 
scriptions of the Brownian dynamics in non-isothermal 
solvents as in the equilibrium case. 

The aim of the present work is twofold. First, to es¬ 
tablish the fluctuating hydrodynamic equations of mo¬ 
tion for a non-isothermal solvent, and secondly to derive 
the coarse-grained description pertaining to a Brownian 
particle suspended therein. Section HIl develops the fluc¬ 
tuating hydrodynamic equations suitable for addressing 
the Brownian motion of a sub-micron sized particle in a 
simple non-isothermal fluid. More precisely, the analysis 
of Sec. U shows that the coupling between momentum, 
temperature and density gives rise to contributions that 
are at most proportional to ei = ATap, where AT is the 
characteristic temperature variation in the system and 
ttp is the isobaric thermal expansion coefficient of the 
solvent. To get a feeling for the numbers involved, con¬ 
sider the paradigmatic example of a hot nano-particle of 
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radius R < 100nm in water [^ . The temperature vari¬ 
ations will usually be bounded by AT < lO^K, so that 
ttp < 10“^ and it is safe to assume ei <C 1. To lead¬ 
ing order, one can thus consider momentum and temper¬ 
ature fluctuations to be independent, and the fluid den¬ 
sity to be constant. Based on these findings, we construct 
the reduced description for the non-equilibrium dynam¬ 
ics of the immersed Brownian particle by eliminating the 
dynamical equations for the solvent fields, in Sec. IIIII 
The particle position X turns out to evolve according 
to a generalized Langevin equation with long-term mem¬ 
ory, whose zero-mean Gaussian noise satisfies a gener¬ 
alized fluctuation-dissipation theorem with a tensorial 
frequency-dependent energy spectrum kBTij{X,u!) that 
implicitly reflects the lack of homogeneity and isotropy in 
the fluid. Finally, in Sec. IlYl we summarize our results, 
leaving a more thorough discussion of the consequences 
on the level of the coarse-grained Langevin dynamics to 
Ref. [13. 


II. FLUCTUATING HYDRODYNAMICS 


The starting point for the following discussion are 
the deterministic hydrodynamic equations describing the 
conservation of mass, momentum and energy in a com¬ 
pressible Newtonian fluid in local thermal equilibrium, 
which occupies the volume V around a suspended Brow¬ 
nian particle of arbitrary shape. Expressing energy in 
terms of the fields T (local temperature) and p (local 
pressure) by means of the local-equilibrium version of 
the first law we have [1^ [13 


dp 


= -pV • V, 


dt 

dt’ „ „ ^ 

p— = -Vp -f V • r - ggz, 


(la) 

(lb) 

(l c) 


where p is the mass density, v the velocity, F the devi- 
atoric stress tensor, Q the heat flux, and Cp the specific 
heat capacity at constant pressure. The gravity force gg 
is directed along the negative 2 ;-axis. We defined the to¬ 
tal derivative -^ = dt+v ■'V to lighten the notation and 
note that, in Eq. (Hi, the temperature variations caused 
by the viscous heating arising from the fluid motion have 
been discarded as second order in the fluxes and there¬ 
fore negligible in comparison with Q. The constitutive 
relations for the deviatoric stress tensor and the heat flux 
read 


r = ?7[Vt;-I-(Vv)’’"]-I-?7b(V • t!)I, (2) 

Q = -kVT, (3) 

where g, rjh, and k are the dynamical shear and bulk vis¬ 
cosities and the heat conductivity, respectively. We also 
introduce the kinematic viscosity v = gj g and the heat 


diffusivity ot = i^l{QCp) (the diffusion coefficients of mo¬ 
mentum and heat), for later convenience. We note that, 
at this stage, all transport coefficients can be thought 
of as spatially varying functions that would have to be 
specified, together with a material law g{p,T), to close 
the system of equations. Having highly incompressible 
solvents like water in mind, we simplify the following dis¬ 
cussion by demanding perfect incompressibility, from the 
outset. Thereby, we forgo the opportunity to faithfully 
discuss very fast processes (faster than the time a sound 
wave needs to travel across a distance of about the par¬ 
ticle size). By moreover neglecting a possible tempera¬ 
ture dependence of the expansion coefficient ttp, which is 
again justifiable for water, the deterministic equations of 
motion for the solvent are closed by the simple material 
relation 


p(r,t) = po 


1 - ap{T{r,t) - To) 


( 4 ) 


where po = p(To) is the density corresponding to the 
reference ambient temperature Tq, and terms of order 
0(ef) were neglected. 

The boundary condition associated to Eqs. m, which 
accounts for the momentum exchange with the suspended 
Brownian particle, is the no-slip condition at the particle 
surface S, i.e. 

v{r,t) = {V{t)) + X r on S (5) 

where (V) and (Jl) are the deterministic translational 
and angular velocity of the particle, respectively [^ . 
The boundary conditions for Eq. m, which describe 
the heat sources maintaining the inhomogeneous temper¬ 
ature field, as well as the particle’s equations of motion, 
are for the moment irrelevant. 

Equations (HD provide the basis for describing the 
deterministic evolution of the coarse-grained non¬ 
equilibrium state of the fluid. Fluctuations about this 
“average state” can be incorporated by adding stochastic 
terms to the stress tensor and the heat flux by substitut¬ 
ing r —>• (5r -I- T and Q ^ 6Q + J, in order to represent 
the random exchange of momentum and energy between 
the hydrodynamic and the omitted microscopic degrees 
of freedom [o,!^- As a consequence, the hydrodynamic 
fields also acquire stochastic contributions according to 
V ^ V + Sv, p ^ p + 6p and T ^ T + ST. Consis¬ 
tency with the local-equilibrium hypothesis sets two con¬ 
straints. First, the probability densities of r and J must 
be Gaussian with mean zero and their variance obeying 
the local fluctuation-dissipation theorem governed by the 
deterministic local values of the temperature field. Sec¬ 
ond, the stochastic equations obeyed by the fluctuating 
fields should be linearized, since non-linear contributions 
to the random fluxes are small by construction. 

We proceed as follows. First, we identify further sub¬ 
leading terms in the deterministic hydrodynamic equa¬ 
tions, which we simplify accordingly. We then insert the 
random contributions to the stress tensor and the heat 
flux to obtain the corresponding stochastic equations of 
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motion. Our analysis closely follows the one used in the 
Rayleigh-Benard problem lea ding to the linearized fluc¬ 
tuating Boussinesq equations [25|. However, an impor¬ 
tant difference concerns the characteristic length scale R 
in the two problems. Namely, in our system the relevant 
length scale is set by the particle size, which we assume 
to be at least several nanometers and less than a micron, 
typically on the order of i? ~ 10“^ m. On this scale, ad- 
vection is much less effective than diffusion in spreading 
momentum and heat in the fluid. The relative magnitude 
of advection and diffusion terms in Eqs. CEi and © is 
estimated as 


O 




~ Re, 


O 


\\aTV^T\J 


~Pe, 


( 6 ) 


where Re = vRjv and Pe = Rev jar are the Reynolds 

I 


and Peclet numbers associated with the particle’s mo¬ 
tion, respectively. To remain consistent with the local 
equilibrium assumption, the characteristic deterministic 
particle velocity (that sets the order of magnitude of the 
fluid velocity v) must remain bounded by the thermal 
velocity Vth yJk^T / {g-pR?'), where R and g-p denote 
the radius and mass density of the particle, respectively. 
In practice, deterministic particle velocities are usually 
much smaller. With pp ~ qq this translates to Re ~ 10“^ 
and Pe ~ 10“^, which implies that advection can be ne¬ 
glected, so that the total derivative d/dt can be replaced 
by the partial derivative dt when acting on v and T, in 
the above equations. 

Under this conditions, and after substituting Eq. o 
into the balance equations (ITal)-(fTB the deterministic 
equations of motion for the solvent degrees of freedom 
become 


OtpOt 


dtT = 


1 - ap{T - To) 


V • V, 


9o 


1 - ap(T - To) 


dtv = -Vp + V - T - go 


QoCp 


l-ap{T-To) dtT = -V -Q 


I - ap{T - To) 
dp 




oipT 


1 - ap{T - To) dt'' 


(7a) 

(7b) 

(7c) 


Notice that, in general, momentum and temperature do 
not evolve independently. In order to understand the 
relative importance of the different terms determining 
such coupling we switch to dimensionless variables: 

- ^ ~R^ ^ ClT 

rR=r, t — = t, V — =v, (8) 

ax R 

TAT = T, f= r, ^OoCpOT^ ^ ^ 

For simplicity we have taken ax and Cp constant, here. 
Equations m in dimensionless form are 


= [l-ei(f-fo)]V-{;, (lOa) 

[I - ei(T - fo)]diV = -V{p + e2z) -h V • f -h eie2^, 

(10b) 

eof dp 


[l-ei{T-To)]diT = -V -Q 


l-ei(f-To) df 

(10c) 


The magnitude of the various terms can be estimated 
by checking the physical values of the dimensionless 
parameters ei = apAT, €2 = R^g/ax, and €3 = 
apQx/(cpRl^), which control the relative magnitude of 
temperature-induced inhomogeneities in density, buoy¬ 
ancy, and pressure-driven heat fluxes, respectively. Tak¬ 
ing R and AT as above, we obtain for water around stan¬ 
dard conditions: ei ~ 10“^, €2 — 10“®, and 63 ~ 10“^. 


This implies that the limit ei,e 2 ,e 3 —>■ 0 reproduces the 
leading behavior of Eqs. while perturbative correc¬ 
tions should follow by expanding the hydrodynamic fields 
in series of these small parameters. To leading order, 
Eqs. (fTHl) then read 

V • -D = 0, (11a) 

dfv = — Vp -I- V • f, (11b) 

dff = -V • g. (11c) 

The condition (llla|) of a divergence-free velocity field 
means that the fluid density can be treated as a con¬ 
stant. Restoring the physical dimensions we obtain our 
final set of deterministic equations 


V • v{r, t) = 0, 


godtv{r,t) 
dtT{r,t) = 


= -V • 


p(r,t)l 


—axV^T{r, t), 


r{r,t) 


(12a) 

(12b) 

(12c) 


where the temperature dependence of the viscosity 
r][T{r,t)] in the deviatoric stress tensor r(r,t) is re¬ 
tained. 

Now we turn to fluctuations and introduce the random 
stress and heat flux into Eqs. (I12bll - (ll2cl) . The resulting 
fluctuating fields obey the equations 
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V • Sv{r, t) = 0, 

Qodt5v{r,t) = -V • 5p{r,t)l - ri[T{r,t)\ {yv{r,t) + VSv{r,tf') + V • T{r,t), 
dtdT{r,t) = -aTV^ST{r,t) - V • J{r,t) . 

Clearly the boundary condition ([S]) becomes 

Sv{r, t) = SV (t) + Sfl{t) X r on 5, 


(13a) 

(13b) 

(13c) 

(14) 


where SV and <5^2 are the stochastic components of the particle velocities. The correlations of r and J are prescribed 
by the local-equilibrium fluctuation-dissipation relations 


{Tij{r,t)Tki{r',t')) = 2r][T{r,t)]kBT{r,t)S{r - r')S{t - t') [SikSp -f SuSjk ), 
{Ji{r,t)Jj{V,t')) = 2aTkBT{r,t)'^6{r - r')5{t - 
{Tij{r,t)Jk{r',t')) = 0 , 


(15a) 

(15b) 

(15c) 


containing only the deterministic part of the temperature field . Consistency with the no-slip boundary condition 
requires that r = 0 on 5 a®. 

Notice that in going from Eq. (I12bl) to Eq. (I13b|) . we have linearized the viscous stress. Indeed, direct insertion of 
fluctuations in Eq. (I12bl) would produce 

r 7 [r + 5r](V('u + (5'u) + (V(u + 5u))T) ~ {p[T] + ^[T]6T){Vv + (Vvf) + r,[T]{VSv + (VSvf) . (16) 


Here we have expanded the viscosity up to first oder in ST 
and dropped the manifestly nonlinear fluctuation term 
proportional to ST'VSv. Equation (I13bl) follows by ne¬ 
glecting the contribution due to temperature variations 
in the viscosity, which is justified by the observation 

f \{dr,/dT)STVv\ \ Ar, , , , 

V \v'^Sv\ p AT Eth ^ ^ 

While the typical relative viscosity variation Arj/rj are of 
the same order as the characteristic relative temperature 
variation AT/Tq, in any coarse-graining volume consist¬ 
ing of N solvent molecules, the spontaneous local tem¬ 
perature fluctuations {ST"^) are small compared to T, 
typically of order 0{T/'/N). Moreover, the boundary 
conditions with the Brownian particle Eqs. ([5]), (fT4l) fix 
the order of magnitude of the deterministic and fluctu¬ 
ating velocity fields to the typical particle velocity Eex 
imposed by external forces and the particle thermal ve¬ 
locity Vth, respectively. As noted above 14x “C Vth is re¬ 
quired to conform with the underlying local equilibrium 
assumption. In typical applications 14x ^ 10~®ms“^ is 
indeed substantially smaller than the thermal velocity 
Eth — 10“^ms“^. 

We thus arrive at the important conclusion that mo¬ 
mentum and temperature are decoupled up to corrections 
of order 0(ei), or smaller. The reason is that only the 
deterministic temperature T{r,t) appears in Eq. (Il3bl) . 
which is fully determined by Eq. (jl2cl) . 

It is interesting, now, to go back to the starting point 
of the present analysis, i.e. the assumption of vanish¬ 


ing solvent compressibility >ct = 0, that permits to as¬ 
sume density variations to arise from temperature hetero¬ 
geneities, alone. A rough estimation of the relative den¬ 
sity variation due to pressure variations can be obtained 
as follows. According to Eqs. (IT^-(fT^. the fluid is 
divergence-free to leading order. The typical magnitude 
of (stationary) velocity and pressure variations can thus 
be approximated, employing the fundamental solution 
of the stationary Stokes equation, by Eh F/rjR and 
Pth ~ F/E?. Eliminating the thermal force F exerted by 
the fluid, we arrive at pth V^th/R — 1 Nm“^. Thus for 
water we get 

o ^ ~ O i^TPth) ^ 10“®, (18) 

which demonstrates the reliability of the assumption 
x'T = 0 on time scales much longer than Ry/Jtrpo — 
10-i°s. 

On the basis of Eqs. m and dni), in the next section 
we derive the generalized Langevin equation for the sus¬ 
pended Brownian particle and its associated noise spec¬ 
trum. 


III. DERIVATION OF THE PARTICLE’S 
GENERALIZED LANGEVIN EQUATION 


We now focus on the non-equilibrium Brownian par¬ 
ticle dynamics. The full state of the fluid-particle sys¬ 
tem is given is terms of the hydrodynamic fields and 
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the particle coordinates, namely the center of mass po¬ 
sition X (t) and the translational and rotational velocity 
Vit) = {Vit)) + SV{t) and n{t) = (f2(t)) -f dn{t), re¬ 
spectively. The latter evolve by Newton’s equations of 
motion 


mV = F^ 

I n = T 


SF^ 

ST- 


F ext, 
■ Text, 


(19a) 

(19b) 


where m is the mass of the particle and I its tensor of 
inertia. The deterministic force and torque exerted by 
the fluid are 

F{t) = — J cr{r,t) ■ n{r) d^r, (20a) 

T{t) = — I r X {cr{r,t) ■ n{r)) d^r, (20b) 

Js 

where n(r) is the inner normal vector field of the particle 
surface S, and 


cr = —pi -b ? 7 r 


( 21 ) 


is the total stress tensor. Analogous definitions hold for 
the random force SF and torque ST, replacing cr by Scr. 
External forces Fext(t) and torques Text(t) may also be 
present. The system of Eqs. (fT^ . (fOl) . (fl^ and (1^ en¬ 
tirely describes the evolution of the fluid and the Brow¬ 
nian particle. Our aim is to eliminate the equations for 
the hydrodynamic fields and reduce Eqs. (HU)-® to a 
generalized Langevin equation for the particle variables 
only. Thus we rewrite Eqs. (HI in the form 


Lb = h + Sh + 


( 22 ) 


where we combine the translational and rotational ve¬ 
locity into the 6-vector b = (V,r2), and we define the 
generalized tensor of inertia 


L = 


and the generalized forces 


m 0 

0 I 


h = 


Sh = 


SF 

ST 


f ext — 


(23) 


By Eqs. (I12bl) . (I12al) . ([5]), and (ITdl) . v{r, t) and p{r, t) are 
linear functionals of b{t') with —oo < t' < t. Therefore, 
in view of Eq. o, the hydrodynamic forces necessarily 
contain a contribution which is a linear functional of b{t') 
with —oo i.e. we can write 


r 

hit) = - / 

J —C 


Z+{X,t-t')- {b{t'))dt'. (24) 


Here Z^{X,t) is a 6 x 6 time-dependent causal friction 
tensor, which depends on the particle position owing to 


the non-constant fluid viscosity. We omit this depen¬ 
dence in the following. The very same reasoning applies 
to Sh, but in addition, since (I13bl) is a non-homogeneous 
equation due to the presence of the random stress r, a 
term ^(t) has to be included in order to account for con¬ 
tributions independent of the particle velocity. Hence Sh 
consists of a friction term and a Langevin noise: 


/ OO 

Z+{X,t-t')-Sb{t')dt' Fi{t). 

-OO 


(25) 


In the subsequent derivation we shall derive the statis¬ 
tics of the Langevin noise ^(t) and relate it to the dis¬ 
sipative term h{t). The linearity of the problem sug¬ 
gests to operate in the frequency space. Given a generic 
function of time git), we denote its Fourier transform by 
giuj) = git)e~^‘^^dt. The complex conjugate of giuj) 

will be denoted by g*(w). 

In Fourier space Newton’s equation reads 


- iu;L ■ biuj) = -Z+iuj) ■ biw) -f |(w) -b /e,,t(w), (26) 

where we used the Fourier transformed Eqs. (HI, (E3 

hiuj) = —Z^iw) ■ (bioj)), (27a) 

Shiuj) = -Z+iu) ■ Sbioj) -b ^iw). (27b) 

Note that the deterministic part of the velocity vector 
(6(w)) is set by the external force and thus can be chosen 
arbitrarily. 

We are now in the position to evaluate the statistics 
of the Langevin noise ^(w). We proceed in three steps. 
First, we derive an expression for the real part of the 
friction tensor defined by 


Zijiuj) = Z^Juj) -b (w). 


(28) 


To evaluate the components of the friction tensor we 
make use of the property Z^ = Z^^, hinging only on the 
symmetry of the stress tensor cry = aji 0. We exploit 
the freedom of choosing the boundary condition ([5]) to 
select velocity vectors whose a-th entry is the only non¬ 
zero one, and denote them by “&i(w) —the superscript 
a will also be appended to the corresponding hydrody¬ 
namic fields. Second, we show that Huj) is a Gaussian 
variable with zero mean. Finally we link the noise corre¬ 
lation tensor (^i(a;)^*(w)) to the friction tensor (051) . 

We wish to find an expression in terms of the solution 
to Eqs. (II2bl) —without formally solving the much more 
involved problem represented by the stochastic equations 
(I13b|) — for the quantity 

Z,,iL0){-bM){^b*iuj)) = Z^piL0){b^iu:)){b*piu:)), (29) 

where the equality holds by virtue of the choice of b(a;). 
In Eq. (051) and in the following we apply the Einstein 
summation convention to latin indices only. Also, we 
suppress the function arguments where there is no risk 
of confusion. Equation (09l) reads 
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# {Zt + zt;){%.){^b*) ^ + (“&.) ^K) 


'*^=® [f^V*) f "aijUjd^r + {^n*) / (r x (“cr • n))^ d^r + (“14) f '’cr* d^r + (“^7^) ( {r x a* ■ n))id?r 
k) s ^ s ^ s ^ s 

= [ ii^V*) + i^n*) X r), dV + [ ((“F) + (“n) x r\ dV 

Js Js 

= ^v*°‘a,jnjd^r + ‘'v/cr-jUj dV 

= / (%*“cry) dV+ [ dj {‘^Vi^a*^) d^r 

Jv JV 

[ / i"v*^Vi - ^v*‘‘Vi) dV 

JV JV 

[ +^^*jdj'^Vi) d^r + iguj [ {-^vl^Vi-^v*°‘Vi) d^r 

Jv Jv 


(30) 


i fTIbl l 


= 2 


f (jiapd^r - 2gujlm ( °‘v*^Vid?r, 

Jv Jv 


where in (1501) we used the divergence theorem and in (ISTl) we defined the generalized dissipation tensor 
(j)af 3 iX,r,uj) = r]{r) {di°‘vj(r,uj)d/v*{r,uj) + d/vj{r,uj)dj°‘v*{r,uj)) = (j)*p^{X,r,uj), 


(31) 


(32) 


where the X-dependence of the hydrodynamic fields is not explicitly displayed. Equation (IHTl) is valid whatever the 
magnitude of “6^ and ^b*, in particular when they are unit vectors. With this choice we have 


Zap{X,uj) = 2 / tpap dJ^r — 2gujlnY / °‘v*^Vi(^r. 

Jv Jv 


(33) 


Since Za^ is real by definition, Eqs. dSSl) and (1321) imply that 4>ap = ‘Ppa- Besides, Eq. (1331) has to be invariant under 
exchange of a and fj owing to the symmetry Z^p = Zpa.. Therefore one can eliminate the second term in Eq. (1331) 
and obtain 


Za/3{X,u}) = 2 [ (j}ai3{X,r,uj)d^r. 

Jv 


(34) 


Then we turn to the random force ^(w): 

t27bl 


= {Sh + Zpbj) (“&,) ^ 6hJ-b,) - -h,5b, 

^ — (“Ei) I SaijUj d^r — j {r x {Scr ■ n))id‘^r + SVi / “cTy rij d^r + 542^ j {r x {°‘a ■ n))i d^r 

Js Js Js Js 

= —J [{°V) + X r) .SaijUj d^r + J [SV + Sfl x r) .°‘aijnj d^r 


la.ira 


L 


"‘ViSa^nj drr + / SvJainUj drr 


= - {"Vi5aij) d^r+ / dj {Svi’^a^j) d^r 

Jv Jv 


tT^ . lTabl 


= — / 5aijdj°‘Vid'^r + / “aijdjSvid'^r — 




"VidjTijd'^r 


ijLJjUUiy. 


""VidjTijd^r 


(35) 

(36) 


= / Tijdj'^Vid^r. 

Jv 

In (|35l) we made use of the divergence theorem, and in dMl) of the property 5aijdj°Vi = ’^UijdjSvi, which is a direct 
consequence of the symmetry of cr. We thus have 


^i( bi) — ^a{ba) 


Tij dj 


:d3r 


(37) 
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which shows that ^ is Gaussian with vanishing mean, being the integral of the deterministic quantity dj°‘Vi times the 
zero-mean Gaussian field t. Hence, its correlation matrix suffices to specify the statistics completely. Using (EH) , we 
determine the noise correlation: 

= [ d^r' f d^rdj‘^v,{r,uj)(Tij{r,uj)T^i{r\uj'))di^vl{r\uj') 

Jv Jv 

= 2k-a5{uj - Uj') [ ri{r)T{r)[di'^Vj{r,oj)d/v*[r,uj) + di'^Vj{r,oj)d/v*{r,ui)) d^r (39) 

Jv 

= 2feBi5(a; - w') ( (j)ap{X,r,w)T{r)d^r. (40) 

Jv 

In (1551) we used the Fourier transform of (|15al) . Setting the magnitude of (6 q,) and {b*^) to one, we finally obtain the 
noise correlation tensor in the form 

(?„(X, a;)Q(X, J)) = k^np{X, uj)Z^p{X, uj)6{uj - J), (41) 

where Tap is the frequency-dependent noise temperature defined by the spatial average of the temperature field T(r) 
performed with the dissipation tensor (f)ap 


Tap{X ,U!) = 


J^(j)ap{X,r,uj)T{r)d^r 

Jv(l^ap{X,r,w)d^r 


Summing up, we have arrived at the generalized Langevin equation for the particle velocity 


L ■ b{t) =-J Z{X, t-t')- b{t')dt' + m + 


(42) 


(43) 


where ^{t) is a Gaussian noise having vanishing mean 
and correlations given by Eqs. (|4T]) . (I42|) . In principle one 
should consider Eq. (IT5a . with the appropriate bound¬ 
ary conditions, in order to determine Eq. (HTl) . If the 
heat sources are independent of the particle motion, e.g. 
if they are placed at the outer boundaries of the fluid, the 
temperature field will generally depend on the instanta¬ 
neous particle position. Nevertheless, in most practical 
cases the thermal conductivities of particle and solvent 
(Kp and K, respectively) will be such that the feedback of 
the particle motion onto the temperature field (which is 
of the maximum relative strength 0 (kp/k — 1) close to 
the particle surface) can be treated as a small correction 
to the overall temperature field in the solvent. If the par¬ 
ticle itself acts as the heat source, like in hot Brownian 
motion [2^ [s^l , the temperature field can be calculated 
once and for all in the particle frame —advection terms 
arising from the change of frame can again be neglected. 
Therefore, Eqs (HU), (HU) and (gSl) alone can be taken 
to entirely describe the particle dynamics. The higher- 
order corrections should remain relatively small, if not 
negligible, for all practical purposes. 


IV. DISCUSSION 

We have analyzed the fluctuating hydrodynamic equa¬ 
tions for a Brownian particle suspended in solvents with 


moderate temperature gradients. The main result ob¬ 
tained in Sec. 1 is that, on the scale relevant for the de¬ 
scription of Brownian motion, it is sufficient to consider 
fluctuations of the solvent hydrodynamic fields around 
their local equilibrium state to linear order. In particu¬ 
lar, the solvent velocity and temperature fields are found 
to evolve independently. This is traced back to two con¬ 
ditions. First, heat and momentum diffusion, rather than 
advection, is the dominant transport mode in the fluid, 
as testified by the small Reynolds and Peclet numbers, 
namely Re ~ 10“^, Pe ~ 10“^. Secondly, commonly 
realized temperature gradients induce negligible relative 
density variations of order ei ~ 10“^. The ensuing 
ineffectiveness of the momentum-temperature coupling 
has a remarkable consequence on the Langevin noise 
Namely, its non-equilibrium energy spectrum is only a re¬ 
sult of the spatial inhomogenity of the random stress ten¬ 
sor T, which is governed by the deterministic temperature 
field T{r,t)- Qualitatively, this could have been antici¬ 
pated by recalling that the most important enhancement 
of fluctuations in non-equilibrium fluids is due to (non¬ 
linear) convective couplings that dominate on long wave¬ 
lengths, beyond the Brownian scale [3l|. Yet, the char¬ 
acteristic energy spectrum k^Tj governing the Brownian 
noise shares with them the distinctive features of non¬ 
equilibrium fluctuations [32|, namely their long-ranged 
nature and their dependence on the mechanical trans¬ 
port properties of the fluid, as encoded in the dissipation 






tensor (pap- 

We conclude that an effective theory, which is linear 
in the particle-fluid momentum exchange and disregards 
stochastic fluctuations in the temperature field, is suit¬ 
able to derive a contracted description for non-isothermal 
Brownian suspensions. The consequences of this result 
are more fully explored in Ref. (27|, where the general¬ 
ized Langevin equation (H51) is derived by direct appli¬ 
cation of linear response theory to the momentum of a 
Brownian particle and to (locally equilibrated) distant 
coarse-grained solvent volume elements. 

Another important conclusion of our above calcula¬ 
tions is that the Langevin noise ^ is characterized by 
Gaussian statistics. Non-Gaussian contributions to the 
fluctuations (around the vanishing mean of^ are pecu¬ 
liar features of out-of-equilibrium systems [33| that can 
appear only if non-linear fluctuations are retained in the 
hydrodynamic equations or if the local-equilibrium as¬ 
sumption for the fluid is violated. Numerical simulations 
[34 . l3^ have shown that the velocity fluctuations of a hot 
Brownian particle clearly exhibit Gaussian distributed 
velocity and (when placed in a harmonic potential) posi¬ 
tion fluctuations, even under extreme heating conditions 
and under narrow confinement. They thus corroborate 
our main result that a local-equilibrium linear fluctuat¬ 
ing hydrodynamic theory with Gaussian noise provides a 
universal basis for deriving equations of motion for non- 
isothermal Brownian dynamics. 

The vanishing average of the Langevin noise implies 
that the important phenomenon of thermophoresis is ab¬ 
sent in our discussion of non-isothermal Brownian dy¬ 
namics. To be precise, the stationary distribution for the 
particle position associated with Eq. (H5)) will in general 
be non-uniform in space, owing to the position-dependent 
noise strength [^ . The corresponding weak kinematic 
thermodiffusion should however not be confused with 
what is commonly referred to as thermophoresis [s^ . 
The notion refers to an effect that cannot be derived 
within any generic hydrodynamic approach, not even if 
density variations induced by temperature gradients are 
considered [l^. It is a result of non-equilibrium molec¬ 


ular interaction forces at the particle-fluid interface and 
therefore remains absent for any hydrodynamic boundary 
condition that does not directly couple the velocity and 
temperature gradient at the particle surface [1^ . Clearly, 
this does not preclude the introduction of thermophoresis 
into the hydrodynamic description by hand, e.g. by im¬ 
posing a hydrodynamic slip velocity at the surface [3^ . 

Finally, it is worth recalling the limitations incurred by 
the incompressibility assumption, kt —t 0, made at the 
outset. While it is certainly adequate to describe the par¬ 
ticle motion averaged over the whole frequency spectrum, 
it is bound to break down at high frequencies, comparable 
to the inverse time needed by a sound wave to propagate 
over a distance comparable to the particle radius. A cor¬ 
responding effect that is well known in the equilibrium 
theory of Brownian motion is the failure of the incom¬ 
pressible theory to recover the mean-squared velocity, 
{V^) = ksTo/m, as predicted by energy-equipartition. 
Instead, it predicts = k^Ti^/M with the renormal¬ 
ized mass M > m accounting for the added inertia of 
the solvent back-flow Q. The discrepancy is straight¬ 
forwardly resolved by observing that the limits t —> 0 
and >ct —>■ 0 do not commute in the evaluation of the 
equal-time velocity autocorrelation function Ac¬ 

cordingly, we expect that the high-frequency predictions 
of our theory will deviate from measurements performed 
with compressible solvents, but presumably only in a fre¬ 
quency range that is currently still difficult to access ex¬ 
perimentally 
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